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In  this  study  the  stability  characteristics  of  shaped-charge  jets  when  exposed  to  axial 
electric  currents  are  investigated.  The  objective  of  this  study  is  to  expand  the  results 
of  previous  analyses  by  Littlefield^®  to  include  high  levels  of  electric  current,  where 
thermal  energy  effects  must  be  included.  Coupling  of  the  magnetohydrodynamic  and 
thermal  characteristics  of  the  flow  is  accomplished  solely  through  the  variation  of 
mechanical,  thermal  and  electrical  properties  with  temperature.  Phase  change  effects 
are  also  included.  The  jet  is  assumed  incompressible  and  perfectly  plastic,  with  the 
Levy-von  Mises  criterion  imposed  to  limit  the  effective  stress.  A  linear  variation  of 
the  effective  stress  with  temperature  is  employed  to  simulate  effects  of  thermal 
softening.  Electrical  resistivity  and  specific  heat  are  permitted  to  vary  linearly  with 
temperature,  coupled  with  associated  jump  values  that  occur  as  the  jet  changes  phase. 
Solutions  to  the  appropriate  base  flow  are  subjected  to  small  axisymmetric  disturbances, 
and  linear  perturbation  theory  is  employed  to  determine  the  time  evolution  of  these 
disturbances.  Perturbations  that  grow  the  fastest  in  magnitude  as  time  progresses  are 
identified  as  the  most  unstable.  Results  of  the  analysis  indicate  that  thermal  effects 
can  dramatically  alter  both  the  base  and  perturbed  liow  fields,  as  well  as  the  growth 
rate  of  perturbations. 


L  INTRODUCTION 

Shaped-charge  jets  are  one  of  the  most  effective  anti-armor  weapons  in  application  today.  Because  of 
their  superior  penetrating  ability,  shaped  charges  are  widely  used  as  warheads  by  all  branches  of  the 
military.  The  penetration  capability  results  from  the  high  velocities  achieved  by  the  jet  tip,  and  from  the 
long  length  facilitated  by  the  axial  velocity  gradient.  However,  penetration  performance  of  the  shaped 
charge  does  have  its  limitations  because  after  a  specified  time  (the  breakup  time),  the  jet  will  particulate 
into  many  fragments  roughly  equal  in  size  and  begin  to  disperse.  The  penetrating  ability  of  the  jet  is 
severely  degraded  after  this  occurs. 

The  particulation  and  breakup  phenomenon  observed  in  shaped  charges  suggested  that  the  stretching 
motion  of  the  jet  may  be  unstable,  and  has  prompted  a  number  of  experimental  and  theoretic^ 
investigations.  Recent  two  dimensional  analyses  by  Curtis',  Pack^  Romero^  and  Littlefield^®  have  been 
successful  in  reproducing  some  of  the  behavior  observed  in  experiments,  indicating  that  the  fragments 
formed  after  particulation  should  have  lengths  roughly  on  the  order  of  die  jet  diameter.  Perturbation 
studies  were  ^so  extrapolated  in  order  to  estimate  jet  breakup  times. 

Recent  theoretical  studies  have  examined  mechanisms  by  which  the  naturally-occurring  instabilities  in 
shaped-charge  jets  might  be  enhanced  in  order  to  promote  premature  breakup.  Once  demonstrated,  these 
mechanisms  would  have  the  potential  for  incorporation  into  armor  packages.  One  mechanism  investigated 
originally  by  Walker®  involves  the  passage  of  an  axial  electric  current  through  the  jet.  Although  this 
technique  was  originally  intended  as  a  means  to  vaporize  the  jet,  it  later  became  apparent  that  axial 
currents  could  be  used  to  promote  the  formation  of  instabilities  'limilar  to  those  seen  in  ideal  plasma 
columns.  This  mechanism  and  others  has  been  investigated  by  Toepfer‘°  and  other  researchers.  Previous 
analyses  by  Littlefield^®  have  demonstrated  that  axial  electric  currents  can  be  used  to  increase  the  growth 
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rate  of  instabilities,  even  when  small  currents  are  employed.  Hie  objective  of  this  study  is  to  expand  the 
results  from  the  previous  research  to  include  high  levels  of  electric  current,  where  thermal  energy  transfer 
mechanisms  become  important. 

A  number  of  assumptions  have  been  made  to  make  the  analysis  tractable  yet  retain  as  realistic  conditions 
as  possible.  The  jet  is  assumed  to  be  incompressible,  infinitely  long,  uniformly  elongating,  and  perfectly 
plastic.  An  axial  electric  current  initially  on  the  surface  of  the  jet  is  applied  at  time  t  =  0,  and  permitted 
to  diffuse  with  time.  This  electric  current  produces  an  azimuthal  magnetic  field,  which  interacts  with  the 
current  to  produce  Lorentz  forces.  The  current  also  results  in  Joule  heating,  increasing  the  internal  energy 
of  the  jet.  Small  disturbances  to  the  idealized  uniform  jet  motion  are  considered,  and  linear  perturbation 
theory  is  employed  to  establish  the  equations  governing  the  time  evolution  of  these  disturbances.  The 
stability  of  the  jet  is  determined  from  solutions  to  these  equations. 

The  arrangement  of  this  paper  may  be  described  approximately  as  follows.  In  Sec.  II,  the  mathematical 
model  and  governing  equations  are  presented.  Solutions  to  these  equations  for  the  idealized  motion  of 
the  jet  are  given  in  Sec.  III.  In  Sec.  IV,  linear  perturbation  theory  is  used  to  develop  the  governing 
equations  for  small  disturbances  to  the  idealized  jet  motion.  In  Sec.  V  solutions  to  the  perturbation 
equations  are  presented  and  discussed.  Finally,  Sec.  VI  contains  some  general  comments  and  conclusions. 

11.  MODEL  AND  GOVERNING  EQUATIONS 

Consider  a  uniformly  extending,  perfectly  plastic  shaped-charge  jet  as  is  shown  in  Fig.  1 .  The  radius  at 
an  arbitrary  point  along  the  jet  boundary  is  denoted  as  rjz,t),  and  the  outward  unit  normal  vector  at  this 
point  as  n.  An  axial  electric  current  l(t),  assumed  to  be  supplied  by  an  external  power  source,  is  applied 
at  time  t  =  0.  In  this  analysis  the  current  is  assumed  constant,  but  it  is  a  straightforward  extension  to 
augment  the  model  for  the  case  of  arbitrary  currents  that  can  vary  with  time. 

Let  the  velocity  be  given  by  V,  its 
components  by  v^,  the  pressure  by 
p,  the  deviatoric  stress  tensor  by 
Ty,  the  current  density  by  J,  the 
magnetic  field  by  B,  the  enthalpy 
by  h,  the  temperature  by  T,  the 
density  by  p,  the  magnetic 
permeability  by  Po-  the  resistivity 
by  X,  the  specific  heat  by  c^,  the 
thermal  conductivity  by  k,  and  the 
yield  strength  by  Y.  Then  the 
governing  equations  in  the  jet  and 
in  the  surrounding  vacuum  are 
given  by 


Fig.  1.  Model  for  Jet  Stability  Calculations. 


pV- v  =  o. 

(1) 

p|^  V  •  VV  j  = -Vp V  •  + J  X  B, 

(2) 

V  .  V/,  J  =  ^-t-V  .  Vp  +X(J  •  +  v  .  (kVT)  +  z,^^, 

(3) 

E  +  VxB  =  Xj,  (4) 

(5) 


V-B  =  0, 

V  X  B  =  PoJ, 


(6) 


and 


V 


(7) 


Equations  ( 1 )  -  (3)  are  the  incompressible  continuity,  momentum  and  energy  equations,  respectively;  Eq. 
(4)  is  Ohm’s  law.  and  Eqs.  (S)  -  (7)  are  Maxwell’s  equations.  Several  other  relationships  are  needed  to 
obtain  closure.  'The  tem^rature  and  enthalpy  are  related  by 


aA_ 

97’"'^'” 


(8) 


and  the  constitutive  model  assumed  for  the  jet  is  a  perfectly  plastic  model,  with  the  effective  stress 
satisfying  the  von  Mises  criterion.  The  deviatoric  stress  depends  only  on  the  velocity  and  the  yield 
strength  under  these  conditions,  and  may  be  expressed  as 


dxj^dxi  J 


(9) 


where  p  is  a  scalar  that  satisfies 

2y2 


XT  =- 

«  «  3 


(10) 


Assuming  the  physical  properties  are  known,  Eqs.  (1)  -  (8)  when  cast  in  cylindrical  coordinates  represent 
a  system  of  10  scalar  equations  for  10  scalar  unknowns:  v„  p,  Ep  7,,  h,  and  T. 


When  large  changes  in  thermal  energy  occur,  several  of  the  physical  properties  that  appear  in 
Eqs.  (1)  -  (10)  can  change  their  values  by  several  orders  of  magnitude.  The  changes  must  be  considered 
in  order  to  make  accurate  predictions  of  the  jet  motion.  In  the  present  study  Y,  k,  X  were  all  pennitted 
to  vary  with  temperature  according  to  the  relationship 


|P,  +  P,7'  T<T„ 

^  lP3+P4r 


(11) 


where  the  Pi’s  are  constants,  T„  is  the  melting  temperature  and /is  th*^  specified  physical  property.  In 
the  phase  change  regime,  properties  were  linearly  interpolated  as  a  function  of  enthalpy  from  their  values 
at  incipient  to  total  melt.  Variations  of  properties  with  resjiect  to  other  state  variables,  such  as  pressure 
or  equivalent  plastic  strain,  were  not  considered.  Although  the  density  certainly  changes  with  temperature, 
its  variation  is  not  as  significant  as  with  the  other  properties  considered  in  Uiis  analysis.  Furthermore, 
the  inclusion  of  variable  density  complicates  the  analysis  considerably  and  will  therefore  be  deferred  to 
future  study. 

Several  of  the  physical  properties  listed  in  Eqs.  (1)  -  (10)  are  discontinuous  at  the  interface  between  the 
jet  boundary  and  the  surrounding  vacuum.  The  density,  for  instance,  changes  suddenly  from  its  value  in 
the  jet  to  zero  in  the  surrounding  vacuum.  The  procedure  for  treating  this  boundary  that  insures  the 
conservation  equations  are  still  satisfied  is  to  integrate  Eqs.  (1)  -  (10)  across  a  layer  of  thickness  5  in  a 
direction  normal  to  the  boundary,  and  then  allow  5  to  approach  zero".  This  procedure  yields  the  correct 
boundary  conditions,  which  are  given  by 


(12) 

-nlp]  +  n  •  IXy]  =  0, 

(13) 

n  •  [kVTl-n  •  [Vp]  +  n  •  [V  •  xJ  =  0, 

(14) 

n  •  [B]  =  0. 


(15) 


and 


n  X  [E]  =  -n  •  v[B], 


(16) 


where  []  denotes  the  change  in  the  specified  quantity  across  the  boundary,  and  v  is  the  velocity  of  the 
boundary.  The  outward  unit  normal  vector  n  may  be  determined  as 


V(r-rt) 

|V(r-r*)r 


(17) 


Equations  (1)  -  (10),  together  with  the  boundary  conditions  given  in  Eqs.  (12)  -  (17),  represent  a  system 
of  coupled  differential  equations  to  calculate  the  time  evolution  of  the  jet  motion.  In  their  most  general 
form,  this  set  of  equations  is  formidable  and  quite  difficult  to  solve.  However,  simplified  solutions  which 
examine  the  stability  of  the  idealized  motion  of  the  jet  may  be  performed  with  a  reduced  set  of  equations. 
Solutions  to  this  reduced  set  of  equations  allow  easier  identification  of  the  physical  principles  that 
contribute  to  the  instability  and  bre^up  of  shaped-charge  jets. 

III.  IDEALIZED  MOTION  OF  THE  JET 

In  the  idealized  motion  of  a  shaped-charge  jet,  the  jet  stretches  uniformly  with  constant  tip  and  tail 
velocities  and  V^n-  The  axial  strain  rate  Ti  is  therefore  T)  =  ( -  V„,,)/L,  where  L  is  the  instantaneous 
length  of  the  jet.  Since  L  increases  as  time  progresses,  t)  decreases.  When  cast  in  a  reference  frame  that 
moves  with  the  center  of  mass,  so  that  v,  =  0  and  z  =  0  at  the  mass  center,  these  kinematic  considerations 
yield  a  velocity  distribution  given  by 

Vrf  =  tl2;  v^--r\rn,  (18) 


where  the  subscript  0  denotes  the  idealized  motion  of  the  jet.  Similarly,  the  length  and  axial  strain  rate 
are  given  by 

=  L=L„x,  (19) 


where  the  subscript  o  denotes  initial  values  and  T  =  1  +  Ti„t.  The  pressure,  enthalpy  and  magnetic  field 
in  the  jet  are  determined  by  solving  appropriate  conservation  equations.  These  variables  may  be  cast  in 
dimensionless  form  by  letting  po  =  fio  =  h(/h„,  To  =  {To-  T^V(T„  -  T„),  So  =  ZTCoBypof .  and 

r  =  r'Jz/a,  where  a  is  the  initial  radius,  I  is  the  electric  current,  is  the  enthalpy  at  melting,  and  B  is  the 
azimuthal  component  to  B.  Use  of  these  dimensionless  variables  in  the  governing  equations  yields 


=  0, 


(20) 


8t 


dpo  EA  .f 


dr  ^  f  j 


2Lii.fr 
pafdfy'^dr  ) 


=  0. 


and 


dSp  So  ^oY 

9x  2x  \  dr  ^  f  ) 


=  0, 


(21) 


(22) 


where  ?o  =  YfJY{T,),  ^  and  ^o  =  Plastic  work  has  been  neglected  in  this  formulation 

because  it  is  .>mall  for  the  problem  at  hand.  The  dimensionless  constants  appearing  in  these  equations 
are 


3p^  Cp^To)[\+i 

'  YiTo)  ’  ^:epr\\a*'  k{TMo  ' 

*'W- 

where  c,  is  the  slope  of  the  specific  heat-temperature  curve  for  the  solid.  Physically,  Cl  represents  a  ratio 
of  inertial  to  plastic  forces,  A  a  ratio  of  electromagnetic  to  inertial  forces,  P  a  ratio  of  "plastic"  to  thermal 
diffusion,  E  a  ratio  of  kinetic  energy  to  sensible  heat,  and  R  a  ratio  of  magnetic  convection  to  diffusion. 
Furthermore,  when  the  physical  properties  are  cast  in  dimensionless  form  using  Eq.  (10),  two  additional 
constants  arise  and  are  given  by 


h„  h. 


where  h,p  is  the  latent  heat  of  fusion  and  /i„  is  the  enthalpy  at  vaporization. 

Boundary  and  initial  conditions  are  required  to  solve  Eqs.  (20)  -  (22).  The  boundary  conditions  are 
determined  from  the  conditions  in  Eqs.  (11)-  (17).  which  yield 


Yq  f  BTq  _  jr  f~ 
SQ=y1', 


Q 


dr 


(23) 


and  are  evaluated  at  f  =  1 .  Initially  the  jet  was  assumed  to  be  isothermal  with  no  electric  current  flowing, 
so  po,  fio>  and  Sq  were  all  set  to  zero  at  T  =  I . 


Shown  in  Fig.  2  is  the  dimensionless 
temperature,  pressure  and  magnetic 
field  distribution  as  a  function  of 
normalized  radius  in  the  jet  at  three 
different  times.  The  parameter  A  was 
chosen  as  100,  whereas  the  remaining 
parameters  were  assigned  values 
typical  for  shaped-charge  jets,  and  are 
listed  in  Table  1 .  While  at  early  times 
significant  surface  heating  occurs, 
central  portions  of  the  jet  begin  to  rise 
in  temperature  at  later  times.  This 
body  heating  is  a  result  of  the 
temperature  variation  in  resistivity, 
which  causes  the  current  density  to 
rise  in  the  center  of  the  jet  where  the 
resistance  is  lower.  The  net  result  is  a 
quasi-uniform  temperature  distribu¬ 
tion  so  that  the  entire  cross  section  of 
the  jet  melts  almost  simultaneously. 


Fig.  2.  Normalized  Temperature,  Pressure  and  Magnetic 
Field  Distribution  vs.  Normalized  Radius. 


In  the  vacuum  surrounding  the  jet,  the  governing  equations 
presented  in  the  previous  section  are  all  identically  satisfied 
except  for  Eqs.  (5)  and  (6).  Furthermore,  since  J  s  0  in  the 
vacuum,  these  equations  reduce  to 


V%  =  0, 

which  yields  the  solution  Sq  =  ^lr. 


Table  1.  Numerical  Values  of 
Dimensionless  Parameters 


£2  0.879 

R  2.13 

E  5.07x10-3 

P  2880 

S  1.50 


Q  4.02 


(24) 


IV.  FIRST-ORDER  EQUATIONS 


Smal]  axisymmetric  disturbances  to  the  uniform  jet  motion  were  then  considered  using  perturbation 
analysis.  Application  of  perturbation  theory  to  the  governing  equations  yields  a  set  of  linear  partial 
differential  equations  that  describe  the  time  evolution  of  disturbances  to  the  jet  motion.  Disturbances 
that  continue  to  grow  as  time  progresses  are  identified  as  unstable.  The  fastest  growing  instability  is  an 
indication  of  the  most  prevalent  disturbance  in  the  actual  jet  motion. 

The  algebraic  details  of  the  derivation  for  the  perturbation  equations  are  complicated  and  quite 
cumbersome,  and  for  brevity  will  not  be  repeated  here.  The  equations  are  formed  by  expanding  the 
dependent  variables  in  a  Taylor  series  fore,  keeping  only  the  terms  are  of  order  e'.  Furthermore,  previous 
analysis  has  demonstrated  that  transformation  of  the  governing  equations  into  the  lagrangian  coordinates 

f  =  rVx/a;  z=zlax  (25) 

permits  the  axial  dependence  of  dependent  variables  to  be  Fourier  analyzed.  Let 


Vw 

= 

Ti„av,exp(ikz“); 

= 

expiikz)-. 

Pi 

= 

iVV  exp(ikz)/4ic^a^; 

hi 

= 

h„fi  expiikz)-. 

= 

pjBexpiikzVlua-, 

= 

{T„-T^)T  expiikz)-. 

Yi 

= 

Yj  expiikz)-. 

= 

Xj^expiikz)-, 

ki 

k„lc  expiikz)-. 

where  k  is  the  axial  wavenumber,  i  =  V-1,  and  the  1  subscript  denotes  a  first  order  quantity.  Then  the 
differential  equations  governing  the  motion  of  small  disturbances  are  given  by 


(27) 
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dx 


dp  J_ 
df  Q 
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In  the  vacuum  surrounding  the  jet.  solutions  to  the  first  order  equations  require  5  =  0  in  the  vacuum.  This 
result  also  follows  directly  from  Ampere’s  law,  since  the  total  current  enclosed  by  the  jet  does  not  change 
when  an  axisymmetric  disturbance  is  applied. 

Several  boundary  and  initial  conditions  are  required  to  solve  Eqs.  (26)  •  (30).  The  initial  conditions  have 
been  presented  elsewhere^  and  are  still  valid  in  the  present  analysis  if  h,  is  assumed  zero  initially. 
Calculations  demonstrated  that  the  initial  value  of  h,  has  a  negligible  effect  on  the  results  anyway,  so  in 
the  present  study  the  initial  v«iiue  of  h/  was  set  to  zero.  The  boundary  conditions  result  from  Eqs.  (11)- 
(17),  and  are  derived  using  a  procedure  similar  to  the  one  used  to  determine  the  governing  perturbation 
equations.  These  conditions  are  given  by 
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and  are  evaluated  at  f  =  1.  In  these  equations,  the  variable  is  the  amplitude  of  the  disturbance  to  the 
radius  of  the  jet,  defined  as 

r^i  =ar^  exp(i7:f)/x''^.  (37) 

The  numerical  value  of  has  special  significance  since  it  is  a  dimensionless  measure  of  the  magnitude 
of  perturbations  to  the  jet  boundary.  Therefore,  disturbances  that  maximize  the  value  of  fj,  are  readily 
identified  as  the  most  unstable  perturbations. 

V.  SOLUTIONS  TO  PERTURBA-nON  EQUA’HONS  AND  STABILITY 
CHARACTERISTICS  OF  THE  JET 

Numerical  solutions  to  the  perturbation  equations  presented  in  Sec.  IV  were  calculated  for  a  wide  range 
of  electric  currents  and  disturbance  wavelengths.  As  a  result  of  the  lagrangian  transformation  applied  to 
the  governing  equations,  the  wavelength  is  time  dependent  and  related  to  the  wavenumber  as  =  Inatlk, 
so  initially  short  wavelengths  become  longer  and  stretch  at  the  same  rate  as  the  jet  length.  Large  electric 
currents  were  applied  in  the  model  by  increasing  the  value  of  the  parameter  A,  in  order  to  in'  estigate  the 
effects  of  large  internal  energy  changes  on  the  jet  stability. 


Shown  in  Fig.  3  is  the  normalized  amplitude  of  the  perturbed  radius  as  a  function  of  normalized  time 
"C  for  several  different  wavelengths.  The  value  of  A  was  1.0,  and  the  values  of  the  other  dimensionless 
parameters  were  as  given  in  Table  1 .  For  this  value  of  A.  therm^  effects  are  uncoupled  from  the  momentum 
and  magnetic  fields  since  the  maximum  value  of  Rq  reaches  only  0.03  when  t  =  1.  Consequently,  the 
values  of  f*  overlay  almost  exactly  the  values  obtained  from  an  isothermal  model®,  where  energy  balance 
effects  were  neglected.  It  is  evident  from  the  figure  that  disturbances  corresponding  to  the  wavenumber 
k  =  K  grow  at  the  largest  rate  for  most  times.  This  wavenumber  corresponds  to  an  initial  disturbance 
wavelength  on  the  order  of  the  jet  diameter.  Furthermore,  disturbances  grow  faster  than  the  case  where 
no  current  is  applied  (A  =  0).  Extrapolation  of  this  result  to  large  amplitude  disturbances  would  indicate 
that  when  low  levels  of  electric  current  are  applied,  jets  should  particulate  into  fragments  with  lengths 
roughly  equal  to  or  slightly  greater  than  the  jet  diameter. 

The  contribution  of  hydrodynamic  and  electromag¬ 
netic  effects  to  instabilities  in  jets  under  isothermal 
conditions  has  been  discussed  at  length  elsewhere®"’, 
and  for  brevity  will  be  only  summarily  discussed  here. 

The  hydrodynamic  instability  results  from  an  axial 
stress  that  opposes  the  stretching  motion  of  the  jet, 
which  is  positive  and  constant  for  the  idealized  jet 
motion.  When  combined  with  the  reduction  in 
cross-sectional  area  in  regions  where  necking  occurs, 
this  constant  axial  stress  results  in  a  net  force  that 
always  tends  to  accelerate  material  out  of  the  neck, 
thus  resulting  in  an  unconditional  instability. 

Moreover,  in  the  absence  of  any  stabilizing  effects, 
this  criterion  would  indicate  that  very  short 
wavelength  disturbances  should  be  the  most  unstable. 

However,  the  short  wavelength  disturbances  are 
partially  stabilized  by  stress  enhancement  that  occurs 
in  the  necking  regions.  The  most  unstable 
hydrodynamic  disturbances,  therefore,  have 
moderate  wavelengths.  The  electromagnetic 
contribution  to  the  instability  occurs  because  the  axial 
current  density  tends  to  increase  in  the  necked 
regions,  resulting  in  an  increase  in  magnetic  pressure. 

This  pressure  increase  also  tends  to  accelerate 
material  out  of  necking  regions,  and  results  in 
unconditional  instability. 

When  the  electric  current  is  increased,  the  behavior  of  instabilities  in  the  jet  change  both  in  magnitude 
and  character.  Shown  in  Fig.  4  is  the  normalized  amplitude  of  the  perturbed  radius  as  a  function  of 
normalized  time  for  disturbances  of  several  different  wavenumbers.  The  value  of  A  was  100,  and  the 
values  of  the  other  dimensionless  parameters  were  as  given  in  Table  1.  For  this  value  of  A,  incipient 
melt  initiates  at  the  surface  at  about  x  =  0.2,  and  complete  melt  begins  to  occur  at  x  =  0.45.  It  is  evident 
from  the  figure  that  larger  currents  result  in  dramatic  increases  in  the  growth  rate  of  instabilities. 
Furthermore,  the  shortest  wavelength  disturbances  appear  to  grow  even  faster  when  A  is  increased.  When 
X  =  0.8,  for  example.  Fig.  4  indicates  that  disturbances  corresponding  to  k  =  271  have  the  largest  perturbed 
radius.  This  wavenumber  corresponds  to  a  wavelength  equal  to  roughly  half  the  jet  diameter.  Moreover, 
it  is  apparent  from  the  slopes  of  f*  that  perturbations  corresponding  tok  =  4n  grow  at  an  even  faster  rate. 
This  seems  to  indicate  that  the  shortest  wavelength  perturbations  will  eventually  become  the  most  unstable 
disturbances. 

The  changes  observed  in  the  growth  rate  of  when  A  is  increased  are  a  result  of  the  additional  importance 
of  electromagnetic  instabilities.  As  the  current  increases,  the  temperature  increases  and  the  jet  material 
softens,  until  the  yield  strength  is  essentially  zero  at  the  melting  point.  Since  hydrodynamic  instabilities 
in  the  jet  are  a  direct  result  of  material  strength,  hydrodynamic  effects  do  not  contribute  significantly  to 
instability  after  incipient  melting.  The  electromagnetic  instabilities,  on  the  other  hand,  become 
increasingly  dominant.  Furthermore,  it  is  evident  that  the  shortest  wavelength  disturbances  should  grow 
faster  when  electromagnetic  instabilities  are  prevalent,  since  the  axial  magnetic  pressure  gradient  that 
activates  the  instability  becomes  larger  as  the  wavelength  diminishes.  This  phenomenon  is  similar  in 
nature  to  the  electromagnetic  instabilities  seen  in  ideal  plasma  columns,  where  the  shortest  wavelength 


Fig.  3.  Normalized  amplitude  f,,  vs. 
normalized  time  for  several  values  of  k.  The 
value  of  A  was  1.0,  and  the  values  of  the  other 
dimensionless  parameters  are  given  in 
Table  1. 


disturbances  possess  the  largest  growth  rates".  Results  in  Figs.  3  and  4  are  in  qualitative  agreement  with 
this  trend,  indicating  that  the  wavelength  of  the  most  unstable  perturbation  decreases  with  increasing 
current. 

It  is  noteworthy  to  compare  the  results  from  the  present  analysis  to  the  results  when  isothermal  conditions 
are  assumed.  Shown  in  Fig.  5  is  f*  as  a  function  of  T  at  several  wavenumbers  calculated  using  an  isothermal 
model**.  The  value  of  A  was  100,  and  the  remaining  parameters  are  given  in  Table  1.  A  comparison  of 
Figs.  4  and  5  indicates  that  the  isothermal  model  underpredicts  the  late  time  growth  rate  of  all  disturbances, 
but  particularly  the  short  wavelength  perturbations.  This  reduction  in  growth  rate  is  a  result  of  the  absence 
of  thermal  softening  in  the  isothermal  model.  Without  thermal  softening,  the  jet  is  partially  stabilized 
by  the  axial  stress  enhancement  that  occurs  in  necked  regions.  Moreover,  this  stress  enhancement  effect 
increases  in  magnitude  as  the  wavelength  diminishes,  which  explains  why  the  short  wavelength 
disturbances  are  stabilized  considerably  when  thermal  softening  is  omitted.  The  magnitude  of  this 
stabilizing  effect  appears  to  be  considerable  and  is  indeed  a  surprising  result,  considering  that  the  plastic 
forces  were  about  two  orders  of  magnitude  smaller  than  the  electromagnetic  forces  in  this  example. 


Fig.  4.  Normalized  amplitude  vs.  Fig.  5.  Normalized  amplitude  vs. 

normalized  time  for  several  values  of  k.  The  normalized  time  under  isothermal 

value  of  A  was  100,  and  the  values  of  the  other  conditions  for  several  values  of  k.  The 

dimensionless  parameters  are  given  in  Table  1.  value  of  A  was  100,  and  the  values  of  the 

other  dimensionless  parameters  are  given 
in  Table  1. 


VI.  CONCLUSIONS 

Large  axial  electric  currents  can  have  a  powerful  disrupting  effect  on  shaped-charge  jets.  Results  from 
the  perturbation  calculations  indicate  that  growth  rates  in  the  perturbed  radius  can  be  increased  by  a  factor 
of  10'*  - 10^  over  their  normal  values  when  A  =  100.  Using  typical  numerical  values  for  strain  rate,  radius 
and  density  in  copper  shaped-charge  jets,  A  =  100  corresponds  to  an  axial  electric  current  of  about  400 
kA.  Although  this  is  a  large  electric  v;urrent,  it  is  definitely  realizable  and  within  the  realm  of  existing 
pulsed-power  technology. 
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